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Abstract: Local metal interconnection is widely used in modern complementary metal oxide semiconductor (CMOS) 
technology. The most frequently used local materials are some heavy metals, such as tungsten (W) or copper (Cu). 
It’s well known that single event upset (SEU) could occur in a CMOS under neutron exposure. In this paper the 
rectangular parallelepiped (RPP) method is used to investigate the SEU response of a typical CMOS. SEU induced 
by 1~14 MeV neutrons are simulated with Geant4 and the cross sections are calculated. The results show that only 
in the structure with W, secondary particle a is created and SEUs are generated when the energy of neutron is less 
than 4 MeV. 
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1 Introduction 

Since some heavy metals, such as copper and tung¬ 
sten, are better conductors than aluminium, metal com¬ 
ponents for interconnection of semiconductor chips can 
be smaller in size with them, and consume less en¬ 
ergy when electricity passes through them while keeping 
higher-performance. Therefore, in semiconductor inte¬ 
grated circuits, tungsten and copper are widely used. 

In complementary metal oxide semiconductor 
(CMOS), incident neutrons, protons, and heavy-ions 
can induce single event upset (SEU) in its sensitive vol¬ 
ume (SV). This paper only focus on neutron induced 
single event upset. In silicon, secondary particles cre¬ 
ated by nuclear interactions of incident neutron can 
range from H to Si, eventually P, including isotopes [1]. 
According to the primary neutron energy and the de¬ 
tails of the interaction, all these particles have variable 
energies, linear energy transfers (LETs) and free paths. 
SEU is caused by the electrical changes associated to the 
electron-hole pairs generated along the ionizing particle 
track. When the number of electron-hole pairs is more 
than the critical charge in the SV of the CMOS, an SEU 
will occur. 

The introduce of heavy metal in CMOS may affect 
the SEU rates of neutron. To study the effect of heavy 
metal quantitatively, we use tungsten interconnection as 
an example. In this paper, a rectangular parallelepiped 
model of a CMOS is built, then the SEU rates induced 
by l~14MeV neutrons are simulated with Geant4 [2]. 


The results of SEU cross sections with and without tung¬ 
sten are compared and the genesis of the difference is 
explored. 

2 Method 

To simplify the structure, the rectangular paral¬ 
lelepiped (RPP) method [3J is chosen to do the sim¬ 
ulation in this work. The method, which simplified 
the structure by using a collection of rectangular paral¬ 
lelepiped, is much easier to model the structure of CMOS 
with Geant4 than a detail description. Deposited energy 
is used to determine an upset instead of total charge in 
the SV because the number of electron-hole pairs is too 
large to simulate. An upset will occur if and only if the 
deposited energy is more than its critical value (E c ), cor¬ 
responding to the critical charge (Q c ). 

Two multi-layer planar targets are used to investi¬ 
gate the SEU response of a typical CMOS with a multi¬ 
layer metallization system, as shown in Fig. |TJ The 
sensitive volume in the silicon for these structures is 
a 3.16 x 3.16 x 2 fim 3 rectangular region which is lo¬ 
cated beneath the metallization stack, where the created 
electron-hole pairs can cause SEU. 

There is only one layer has difference between the 
two structures. In one structure, this layer is composed 
of tungsten, as shown in (a) of Fig. [TJ which is commonly 
used in integrated circuits to provide electrical connec¬ 
tions between layers of metallization or in contact with 
the underlying silicon j3], while in the other structure it 
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is silicon dioxide, referred to (b) of the figure. This layer 
is approximately 1.4 fim above the sensitive volume with 


0.6 p,m in thickness. Characters of other layers are also 
shown in Fig. [I] Effects of the tungsten layer on the SEU 
cross section are presented in this paper. 




(a) 


Ti 0.1pm 


] TiN 0.1pm 


(b) 


Fig. 1. Structures of the scaled CMOS without (a) and with (b) a tungsten layer. Only one layer has difference 
between the two plots. 


Energy deposition in the SV of the two structures 
is simulated with different mono-energetic beam of neu¬ 
tron. Total 6 x 10 8 neutrons are generated and perpen¬ 
dicularly inject into the top of the CMOS. The elastic 
scattering, inelastic scattering, neutron capture process 
and ionization process for the secondary ions are taken 
into account in our code. 

In silicon, an electron-hole pair creation normally 
needs around 3.6 eV deposited energy [5], so the rela¬ 
tionship between E c and Q c is 


3 Results 


1~14 MeV neutrons are simulated with Geant4 in 
this paper and the cross sections are calculated, respec¬ 
tively. The simulation results of the two structures are 
shown in Fig. [2] 


E c (MeV) = 22.5 xQ c (pC). 

The critical energy of our CMOS is set as 0.69 MeV [6]. 

The cross section of single event upset can be ob¬ 
tained with the following relation: 

7—i \ n SE u 
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where E„ is the energy of incident neutron, Useu is the 
number of SEU events with E dep > E c , E dep is deposited 
energy in the SV. And (f> is the fluence of the beam, which 
is calculated with 



where N is the total number of incident neutron simu¬ 
lated and A is the irradiated area. It’s worth to note 
that the cross section depends on E n and E c , which is 
target geometry dependence. 
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(b) 

Fig. 2. SEU cross sections of the two structures 
simulated with different energy neutron. The blue 
filled circles denote the cross section of the struc¬ 
ture with W, while the red filled rectangles for the 
one without W. The data in (b) are as same as 
(a), except shown with log-Y axis. The points are 
shifted slightly from their values to avoid overlap¬ 
ping. 


As shown in Fig. [2] the SEU cross sections of the 
two structures are almost the same when the energy of 
neutron is higher than 4 MeV. In other words, the tung¬ 
sten layer almost has no effect on the SEU rates for high 
energy neutron (> 4 MeV). But when the energy of neu¬ 
tron is less than 4 MeV, no SEU has been seen in the 
structure without a tungsten layer, indicating the tung¬ 
sten layer enhanced the SEU probability in low energy 
area of neutron (< 4 MeV). 


1 MeV neutron in the structures with and without W are 
checked and shown in Fig. [3] In the structure without a 
tungsten layer, no event has been seen with the deposited 
energy in the sensitive volume greater than the E c , while 
there are 3 out of 6 x 10 s events with deposited energies 
greater than the E c value in our simulation. The domi¬ 
nated contribution of energy deposition in the SV is from 
secondary particles created by neutrons through nuclear 
interactions. The total deposited energy of each kind of 
secondary particle in the SV from the simulation result 
is obtained and filled in a histogram as shown in Fig. [4j 
which is normalized into one neutron. The difference be¬ 
tween the two series of data indicates that the secondary 
particles in the two structures are different because of 
the existence of the tungsten layer. 
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Fig. 3. The distribution of deposited energies in 
the SV. The blue filled circles denote the results 
for the structure with W, while the red filled rect¬ 
angles for the one without W. Total 6x10 s 1-MeV 
neutrons are simulated. The vertical red line is for 
E c = 0.69 MeV. 

To explore the increase of SEU rates in the struc¬ 
ture with W, the energy deposition distributions of 


Fig. 4. Species of secondary particles and the cor¬ 
responding energy deposition normalized into one 
neutron in the SV. The blue filled circles denote 
the results for the structure with W, while the 
red filled rectangles for the one with Si02. Total 

6 x 10 8 1-MeV neutrons are simulated. No a and 

7 has been seen in SV of the structure with Si02. 


Figure [4] shows that secondary particles produced in 
both structures are 28 Si, 29 Si, 30 Si, 16 0, electron, and 
proton. The particle which produced the maximum to¬ 
tal deposited energy is 28 Si in both structures. The dis¬ 
tribution of deposited energy of 28 Si is shown in Fig. [5j 
and that of 29 Si, 30 Si and 16 0 are shown in Fig. [6] Plots 
for deposited energies from electrons, protons and js are 
ignored in this work since they are too few. 
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Fig. 5. Deposited energies of 28 Si in structures 
with and without tungsten layer. Total 6 x 10 8 
1-MeV neutrons are simulated. The blue filled 
circles denote the deposited energy for the struc¬ 
ture with W, while the red filled rectangles for 
the one without W. Detailed description can be 
referred in the text. 



Fig. 6. Deposited energies of 29 Si, 30 Si and ie O in 
structures with and without tungsten layer. Total 
6 x 10 8 1-MeV neutrons are simulated. The red 
filled and open circles denote the deposited energy 
for 29 Si in structures with and without W, while 
the blue filled and open rectangulars for 30 Si with 
and without W, and black filled and open stars 
for ls O with and without W. 


Figure [5] shows that in both structures, a large num¬ 
ber of 28 Si are generated, but the energies deposited by 
them of each event are much less than the critical energy. 
The W layer is replaced by Si0 2 in the other structure, 
so more 28 Si are created in it, which result in more de¬ 
posited energy from 28 Si in the structure with Si0 2 than 
that in the structure with W. The deposited energies of 
the rest particles for each event are also far less than the 
E c in both structures, as shown in Fig. [G] According to 
these figures, the corresponding particles are not the rea¬ 
son of the increase of SEU cross section in the structure 
with W layer. 


Figure [3] shows that the clear difference between the 
secondary particles in two structures is that many as are 
created in the structure with W layer while no a in the 
structure without it. The energy deposited by a of each 
event in the SV of the structure with W layer is shown 
in Fig. 0 



Fig. 7. Deposited energy of alpha in the SV of 
the structure with W layer. Total 6 x 10 8 1-MeV 
neutrons are simulated. The vertical red line is 
for E c = 0.69 MeV. 


Although there are only a small number of a created 
in the sensitive volume, most of them have much larger 
deposited energies than that of the other particles in the 
structure with W, which indicates that events with a 
created have more probability to deposit high energies 
and generate SEUs. 

Extended studies show that when the energy of neu¬ 
tron is higher than 4 MeV, no a is seen in the struc¬ 
ture without W. But more kinds of secondary particles 
are created and the probability to deposit enough energy 
causing SEUs increases for higher energy neutron. The 
SEU cross sections of the two structures are similar in 
the relatively high energy region, as shown in Fig. [2} 

4 Conclusion 

Single event upset rates induced by 1~14 MeV neu¬ 
trons in a multi-layer metallization CMOS are simulated 
with Geant4 in this work and the cross sections of SEU 
are calculated. No SEU has been seen in the struc¬ 
ture without W when the energy of neutron is less than 
4 MeV, while the cross sections are similar for the struc¬ 
tures with and without W because the secondary parti¬ 
cles other than a have chance to deposit enough energy 
to cause SEUs when the incident neutron energy is higher 
than 4 MeV. 

Detailed studies with 1 MeV neutron indicate that 
the secondary particle a is created with the existence 
of W layer, while no a has been seen in the structure 
without W. Only events with a created can have chance 
to deposit enough energy to cause SEUs with such low 
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